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ABSTRACT 

A primordial cosmic microwave background B-mode is widely considered a "smoking 
gun" signature of an early period of inflationary expansion. However, competing the- 
ories of the origin of structure, including string gases and bouncing cosmologies, also 
produce primordial tensor perturbations that give rise to a B-mode. These models 
can be differentiated by the scale dependence of their tensor spectra: inflation pre- 
dicts a red tilt (tit < 0), string gases and loop quantum cosmology predict a blue tilt 
{tit > 0), while a nonsingular matter bounce gives zero tilt (ny = 0). We perform 
a Bayesian analysis to determine how far \tit\ must deviate from zero before a tilt 
can be detected with current and future B-mode experiments. We find that Planck 
in conjunction with QUIET (II) will decisively detect tit ^ if \tit\ > 0.3, too large 
to distinguish either single field inflation or string gases from the case tit = 0. While 
a future mission like CMBPol will offer improvement, only an ideal satellite mission 
will be capable of providing sufficient Bayesian evidence to distinguish between each 
model considered. 
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1 INTRODUCTION 

Precision cosmology is on the verge of answering deep 
questions about the early Universe. Are the temperature 
fluctuations of the cosmic microwave background (CMB) 
Gaussian? Will we discover primordial tensor modes in 
the CMB? The existence of a large-scale spectrum of ten- 
, sor perturbations is widely considered to be indicative of 
an early period of inflation - a "smoking gun" signa- 
ture of quasi-de S itter expansion ; howe ver, inflation is not 
the only source (|Brandenbergerl l201ll ). Alternative theo- 
ries of the origin of structure, like string gas cosmology 
or bouncing Universe models, are capable of generating 
B-modes, as are topologi cal defects (jSeliak fe Slosarl [20061 : 

iPogosian fc Wv man 200 8]) and bubb le collisions from phase 
transitions (|Jones-Smith et al . 2008) occurring after the big 
bang. These latter two sources have the potential to obscure 
or confuse the signal from inflation or its alternatives, but fu- 
ture observations should be able to dis entangle the sources if 
the primordial signal is large enough dUrrestilla et ak 2008; 

iBaumann fc Zaldarriagal 120091 ; iGarcia-Bellido et al.ll201lT ). 
The purpose of the present paper is to investigate how well 
the different primordial sources of tensor modes - inflation, 
string gases, or bounces - can be resolved and supported. 
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The tensor perturbation spectrum is typically modeled 
as a power law in comoving wavenumber, k, 



P,Ak) = />,.(/.-„)( y 



(1) 



where tit is the tensor spectral index. In slow roll infla- 
tion the amplitude of tensors is determined by the Hub- 
ble parameter, Pu oc H 2 , which, together with H < 
in the Friedmann Universe, furnishes the hallmark predic- 
tion of a red spectrum: tit < 0. Meanwhile, the spectrum 
of metric fluctuations that arises from a string gas in the 
early Universe is characterized by tit > in the string 
frame (|Brandenberger et al.l [2007), due to the growth of 
anisotropic stress as the radiation dominated phase is ap- 
proached. A blue spectrum can also be generated du ring a 
phase of so-called superinflation (|Piao fc Zhand [20041). mo- 
tivat ed, for example, by loop quantum cosmology ( Boiowaldl 
2002). Yet a different picture emerges from fluctuations pro- 
duced during the contracting phase of a bouncing cosmol- 
ogy: a matter dominated contraction with a regular bounce 
gives a scale invari ant spectrum i n the expanding radiation 
domi nated phase (jWandsl 1 19991 ; iFinelli fc Brandenbergerl 
120081 ) . Evidently, the sign of tit is a more likely indicator of 
inflation than the mere presence of tensors. But how large 
must this tilt be in order for current and future probes to 
detect it with confidence? 

The basic question that we seek to answer in this work 
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Figure 1. A simulated 2cr detection of tit 7^ obtained from 
mock Planck data with a fiducial tensor spectrum r = 0.1 and 
n T = 0.3. 



is how far tit must deviate from zero before one can confi- 
dently conclude tit 7^ 0. From the point of view of sam- 
pling statistics, this is a straightforward hypothesis test: 
how large must tit be in order to reject the null hypoth- 
esis, Ho : tit ~ 0? Having analyzed data, d, one determines 
whether the best fit value, ut, is a sufficient number of stan- 
dard deviations, o~ nT u, away from tit — to satisfy the 
chosen significance level, a. Specifically, assuming that tit 
is Gaussian distributed, the p-value must be calculated, 



-/2 



dz ^ 



(2) 



where the test statistic s — \ut\/ '<Jn T \d is the number of 
standard deviations that tit lies away from tit = 0. For 
example, the significance level a = 0.05 corresponds to s = 
1.96, with the conclusion that one risks a 5% chance of falsely 
rejecting the null hypothesis, tit = 0. 

If one were to proceed along these lines, it would be 
determined that the Planck SurveyoiQ will detect ut 7^ 
at the 95% confidence level if \tit\ 0.3 for a tensor/scalar 
ratio r = 0.1 (see Fig. l.jf] What this approach fails to de- 
termine, however, is whether there is a need for allowing 
ut to vary in the first place: is the fit to the data suffi- 
ciently improved to warrant the inclusion of ut as a free 
parameter? The accepted approach to this problem is to ap- 
ply Bayesian inference to the space of competing models. 
As a definition of model, we should have in mind a collec- 



1 http://www.esa.int/planck 

2 This result is obtained without the assumption of Gaussianity 
leading to Eq. p), but instead follows from a Bayesian parameter 
estimation performed on a simulated Planck-precision data set 
with fiducial tensor spectrum r = 0.1 and tlt = 0.3. We will 
discuss our analytical method in detail in the following sections. 



tion of parameters, = (6\, • ■ • , 9 n ), with associated prior 
probabilities tt(#i), • ■ •, ir{9 n ). Here we have two models: Ho 
corresponds to the null hypothesis in which tit = 0, and Hi 
contains ut as a parameter free to vary across some prior 
range, [n<r,mm, nr^ax]. 

Bayesian analysis gives pr eference to the l east complex 
model that best fit s the data dMacKavll2003l: lTrottall2005l: 
Liddle et all |2006| ; iMukheriee fc ParkinsonTbOlOl : iTrottal 



2008). Models with many parameters free to vary across 



large prior ranges, A9j ~ l/n(9j), are considered more 
complex than models with fewer parameters and more re- 
strictive priors. Bayesian selection prefers models that are 
predictive: the prior and posterior parameter widths should 
be comparable. Frequentist significance tests like Eq. ((2| 
do not incorporate this essential aspect of model selection, 
and can sometimes give c onclusions in striking disagreement 
with Bayesian inference (|Lindlevlll95"7l : IShaferlll982l : ITrottal 
2005). This disagreement, known as Lindley's paradox, tends 
to be most prevalen t for detections in the range of two to 
four a (|Trottai r2005V the same threshold at which the sam- 
pling statistics hypothesis test Eq. © will detect ut = 0.3 
with Planck. In this work, we do not solely examine how 
well future probes might constrain tit, an analysis carried 
out in detail by othe rs (| Verde et alj [2006: Zh ao fc Zhand 
120091 : iMa et alJ lioioh : we are instead interested in going 
one step further, and determining whether or not, given 
these constraints, we are warranted in including ut as a 
free parameter. In the next section, we develop the frame- 
work for performing Bayesian model comparison forecasts 
for future experiments , follow ing closely the approac h taken 
by l|Pahud et al.l I2006T ) and (jVardanvan et al.l T2006) to in- 
vestigate the tilt of the scalar spectrum and the curvature 
of the Universe, respectively. 



2 BAYESIAN EVIDENCE 

Methods of Bayesian inference can be applied to the problem 
of model selection in a manner fully consistent with proba- 
bility theory. Given a model, Hi, and a dataset, d, Bayes' 
Theorem gives 

p(d\Hi)n{m) 



p(H t \d) = 



7T(d) 



(3) 



where the quantity p(Hi\d) is the posterior probability of 
the model Hi, and p(d\Hi) is the evidence for model Hi. 
One typically assigns equal prior probabilities ir(Hi) to the 
alternative models, and 7r(d) is independent of the model. 
If we again apply Bayes' theorem, 

p(d\e,Hi)ir(e\Hi) 



P (e\d,Hi) = 



(4) 



p(d\Hi) 

we obtain the posterior probability of the parameters given 
d. The term p(d\8, Hi) is a function of both the parameter 
values and the data; for a given dataset, it is commonly 
called the likelihood of the parameters, C(0\d, Hi). The ev- 
idence from Eq. @ is nothing more than a normalization 
constant for Eq. Q, 



p(d\Hi)= £(O\d,Hi)n{0\Hi)d0. 



(5) 



This expression defines the evidence to be the prior-weighted 
average of the likelihood over the parameter space. 



Tensor Tilt from Primordial B-modes 3 



For the sake of illustration, assume that the parameter 
6j is well constrained by some data, d, so that the posterior 
probability p(9j\d, Hi) is well-peaked within the prior range, 

1 



\Hi) — [H(9j — 9j >U iin)H (#j,max — 9j)\ 



A8i 



(6) 



where H is a step function that enforces the prior range and 
we have chosen a, flat prior on 9j, n(9j) = 1/A9j, within this 
range. Since the posterior distribution is well-peaked, we can 
approximate the integral Eq. (JS| using Laplace's method 
with Eq. (|4]): we simply multiply the height of the un- 
normalized posterior (the numerator in Eq. ©), p(9j jd, Hi), 
by its width, a g , d , 



p(d\m) 



\d,Hi)de 3 



p(6j\d,Hi) x a B . 



(7) 



where 9j denotes the point of maximum likelihood. In terms 
of the likelihood function this becomes, 



p(d\m) 



cie^H^ie^Hi)^, 
c[0:d.H:"' ' 



A0, 



(8) 



While only an approximation, this expression nicely reveals 
the essential ingredients of Bayesian model selection: a high- 
valued maximum likelihood clearly increases the evidence in 
favor of the model, while the Occam factor, /3 — a g .\ d j A9j ^ 
1, penalizes overly complex or poorly predictive models. 
Models with a prior volume much larger than the posterior 
volume, /) < 1, are not considered predictive because they 
can accommodate a wide range of parameter values before 
the data is collected. Complex models with loose priors and 
many free parameters that are well-constrained by the data 
are therefore penalized by the Occam factor, and will conse- 
quently have lower evidence than a simpler, more predictive 
model that fits the data equally well. 

We are now ready to do model selection: we simply com- 
pute the Bayesian evidence for each model and compare. 
Given two competing models, Ho and Hi, this can be done 
via the Bayes factor, 

p(d|tto) 



Bo 



(9) 



p(d|«i) 

A rubric for scoring the signific ance of a mod el is given 
by the well-known Jeffreys' scale ( Jcffr evslll96ll ). The scale 
rates: | ln_Boi| < 1 (indecisive), 1 < | ln_Boi| < 2.5 (substan- 
tial), 2.5 < | In Boi| < 5 (strong), and |lnB i| > 5 (deci- 
sive), with In Boi > (In Boi < 0) favoring Ho {Hi). In this 
work, we will quote results for strong and decisive evidence, 
corresponding to odds ratios of 12:1 and 150:1, respectively. 

While the Laplace method is instructive, it is not suf- 
ficient for an accurate determination of the evidence. While 
evaluating the integral in Eq. ([5]) is computationally de- 
manding, various methods have been applied to problems 
of mo del selection i n astrophysics, includin g nested sam- 
pling |Skilling| [iobl : iMukheriee et al l |200rl) and thermo- 
dyna mic integration (|Q'Ruanaidh et al.ll 19961 : iBeltran et al.l 
120051). Here w e make use of the Savage-Dickey density ratio 
( Dickevl[l97ll ). an exact analytical expression for Boi that 



can be applied whenever the models to be compared are 
nested. Models Ho and Hi share the same cosmological pa- 
rameters, i/>, except for nr which is set to zero in Ho- For 



separable priors, ir(tp\nT, Hi)\n T =o = K(tp\Ho), the Bayes 
factor takes the form 



Boi 



p(n T \d,Hi) 
■k(71t\Hi) 



(10) 



where p(riT\d, Hi)\n T =o is the marginalized posterior prob- 
ability of tit under model Hi, 



p(n T \d,Hi)= £({*/>, n T }\d, Hi) di/>, 



(11) 



evaluated at tit = 0. This quantity can be obtained rel- 
atively easily using Markov Chain Monte Carlo (MCMC) 
techniques, and, in principle, gives Boi as a function of nr- 
This is the function that we seek to determine in this analy- 
sis, for a variety of current and proposed CMB experiments. 
For each experiment, we will obtain projections by generat- 
ing simulated CMB data across a range of values of nr, and 
determine how the evidence for Hi builds as |tit| grows. 



3 CMB SPECTRA FROM FUTURE PROBES 

Primordial perturbations impart inhomogeneities in the 
photon temperature at decoupling, measured today as di- 
rectional anisotropies on the last scattering sphere, 



5T(n) 
T 



= y^,aJmYem(n), 



(12) 



where the multipole moments, a£ m , are complex Gaussian 
random variables with variance {aj^aj, m ,) — Cf T 5 U iS mm i 
in the direction n. Quadrupolar temperature anisotropies 
at decoupling (and again at reionization) are projected into 
anisotropies in the polarization of the CMB an d can be sim- 
ilarly decomposed ( Kamionkowsk i~et al. I [19971 ). 



Tab \- r 



m 1 (£m)ab 



(n 



i B 

t a e „ 



Y, 



(£m)ai)( n )j ) 



(13) 



where the Y, £ ^. ab are electric-type (curl-free) and magnetic- 
type (divergence-free) tensor spherical harmonics, respec- 
tively. The polarization anisotropies are described by the 
correlations, 



l E i 2 \ _ n EE 



|2\ _ r BB 

l t.m. \ I — y->t. , 



and one nonzero cross correlation, 

I T* E \ _ n TE 
\ a lm a lm) — W 



(14) 



(15) 



Primordial density perturbations can be constrained 
by measurements of the temperature and E-mode po- 
larization anisotropies, while primordial gravitational 
waves additionally create a B-mode polarization pa ttern 
l|Kamionkowski et al.lll997l ; ISeliak fc Zaldarriagalll997r i. The 
B-mode signal is therefore a key indicator of primordial ten- 
sors. 

Our projections are based on simulated datasets. Since 
we do not have access to the true distribution of the ai m 's, 
an estimator is formed from their measured values, 



C 



= E 



\ a fm a ^lm\ 
21+1 



(16) 
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where XY = TT, EE, BB, and TE. As the sum of the 
squares of Gaussian random variables, the Cf Y are xt' 
distributed with v = 21 + 1 degrees of freedom. We gen- 
erate simulated data by dr awing the C f Y 's from a X2t+i 
distribution with variances (Knox 1995), 



(AC? 
(AG 



FE\2 _ 



(2£+l)/ sky 
2 

(2£+l)f sky 



Cf x + N* x ) 



(17) 



a 



teV 



+ 



(cf T +Nl T ] 



(cf E + Ni E 



)] 



(18) 



where / s k y is the fraction of sky covered, the Cf Y are the 
theoretical signal spectra, and the N* Y are instrumental 
noise spectra for X = T, E, and B. For experiments with 
multiple frequency channels, c, the full noise spectrum is the 
inverse sum of the individual spectra over the channels, 



(19) 



Assuming a Gaussian beam, the noise spectrum character- 
izes the combined effects of the instrument beam smearing 
and Gaussian white pixel noise, 



at; 



(opL#fwhm) 2 exp 



9 2 

'fwhm 



8 In 2 



(20) 



where crjjf x is the noise per pixel (with cr£ ix = y/2ap ix ), f?f w hm 
is the full width at half maximum (FWHM) of the Gaussian 
beam, and pixel noise from temperature and polarization 
maps is uncorrelated and taken to vanish. 

We obtain projections for the Planck Surveyor, both 
alone and in combination with the ground-based telescopes 
PolarBeaJf] and the Q/U Imaging Experi ment (QUIETfl 
the p roposed CMBPol satellite mission (|Baumann et al.l 
120091 ) . and an ideal satellite. The satellite missions are sensi- 
tive to both the temperature and polarization spectra, while 
the ground-based experiments will measure only polariza- 
tion. The two platforms are also sensitive to B-modes in 
different ranges: satellites are most sensitive to the reion- 
ization hump at £ < 10, while ground-based detectors will 
perform best at intermediate scales, I = 30 — 500, probing 
the anisotropy at decoupling. In Figure 2 we show the S/N 
ratios of the experiments considered in this analysis as a 
function of multipole number, £, for a B-mode signal with 
r = 0.1. 

For Planck, we include three channels with frequencies 
(100 GHz, 143 GHz, 217 GHz) and noise levels per beam 
(<^pix) 2 = (46.25 fiK 2 , 36 fiK 2 , 17.6 /iK 2 ). The F WHM of 
the three channels are ffwhm = (9.5', 7.1', 5.0') (|Planckl 
2006). For PolarBear we consider frequency channels (90 
GHz, 150 GHz, 220 GHz) with (crjL) 2 = (2.6 fiK 2 , 5.8 
[iK 2 , 744 ^K 2 ) and resolutions fwhm = (6.7', 4.0', 2.7'), 
and for QUIET (II) we use frequencies (60 GHz, 90 GHz) 
with (a^) 2 = (0.04 ^K 2 , 0.08 ^K 2 ) and (9 fwhm = (23', 10') 
ijSamtleben et al1l2007l V We assume integration times of 0.5 
and 3 years for these two experiments, respectively. We com- 
bine three channels for CMBPol with frequencies (100 GHz, 



3 http:/ /bolo.berkelcy.cdu/polarbear 

4 http://quiet.uchicago.edu 
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Figure 2. Signal-to-noise ratios of experiments considered in this 
analysis as a function of multipole number, i, based on a fiducial 
model with r = 0.1. 



150 GHz, 220 GHz) and noise levels (q^ Y ) 2 =(729 nK 2 
676 n K 2 , 1600 nK 2 ) and 8 {whm = (8', 5', 3.5') l|Fraisse et all 
2008). In Figure 2, for reference, we also include the S/N ra- 
tio with lensing as the only noise source (green dashed). 
This contamination, arising from the gravitational lens- 
ing by large scale structure of pri mordial E-modes into B- 
modes dZaldarriaga fc Se liak 1998), dominates the signal at 
£ > 100; probes that seek to measure the B-mode signal 
in this range must successfully subtract the lensing contri- 
bution. This effect helps motivate our selection of exper- 
iments: PolarBear's sensitivity lies just below the lensing 
signal, while the full capability of the upgraded QUIET (II) 
and CMB Pol will be more sensitive and require succe ssful 
delensing (|Hu fc Okamotdliool Iffirata fc Seliakll2003l ). We 
assume for these experi ments that lensing is reduced to a 
level of 5 x 10" 8 yK 2 (|Seliak fc Hiratall2004h . Lastly, our 
ideal experiment will suffer from zero instrumental noise or 
beam effects, but will still be limited by cosmic variance, 
(ACi) 2 = 2Cf/{2£ + 1), and the reduced lensing noise. We 
assume sky coverages of / s k y = 0.65, 0.012, 0.04, 0.8, and 
1.0, respectively, for the experiments just listed. We next 
present our analysis and results. 



4 EVIDENCE FOR A TENSOR TILT 

The determination of how Boi varies with tit requires an 
evaluation of p(nT\d,Hi) (c.f. Eq. [9| across a range of nx, 
although, in general, this distribution also depends implic- 
itly on r (a„ T is an increasing function of r.) We first con- 
front the simpler task of obtaining results for Boi with the 
fiducial tensor amplitude fixed at r = 0.1, near the uppe r 
95% C.L. set by WMAP+BAO+i/o (|Komatsu et all 1201 ll ). 
Though limited, this provides an optimistic projection and 
filters out those experiments that fail to provide strong evi- 
dence for tit 7^ even under the most favorable conditions; 
those experiments that do provide strong evidence are ana- 
lyzed further for the case of variable r. 



4.1 Optimistic Projection 

For each experiment listed in the previous section, we obtain 
p(jit\&, Hi) by analyzing simulated data generated from the 
model Hi. We assign the base cosmological parameters the 
fiducial values Q. b h 2 = 0.022, Q c h 2 = 0.105, 6 S = 1.04, 
r = 0.09, A 3 = 2.23 x 10" 9 , n a = 0.97, r = 0.1, and 
simulate a different dataset as tit is incremented in steps 
of \S nT \ — 0.1 across the prior range tit G [— 0.5, 0.5]j. 
We use MCMC to constrain r and tit for each dataset; 
since the base parameters Q.bh 2 , Q, c h 2 , 6 S , r, A s , and n s 
have little effect on the B-mode signal, we only vary r and 
tit within the chains. The theoretical temperature and po- 
larizati on Ci-spectra are generated out to I = 2000 with 
CAM^f] (|Lewis et al.ll2000l ), and the parameter est imation is 
performed using CosmoMQj (|Lewis fc Bridle! 120021 ') . For each 
dataset, convergence is measured across four chains using 
the Gelman-Rubin R statistic. The constraints on r and ut 
are in general correlated, but the parameter uncertainties 
can be minimized by choosing the pivot scale corresponding 
to the multipole at which they become uncorrelated, fc* ~ 



10" 4 Mpc" 



. This pivot scale depends on the data, and will 



be different for each experiment: £* = 10, 30, 65, 150, and 
300 for Planck, Planck+PolarBea r Planck+QUIET (II) 
CMBPol, and the ide al experiment l|Zhao fc Baskaranll2009l ; 
I Zhao fc Zh ang 2009). By constraining r(fc„), we minimize 
the effects of correlations in our comparisons of constraints 
across experiments. 

We present results in Figure 3 and Table 1. In Figure 
3, the data points represent the actual evaluations of lnSoi; 
the curves are obtained via quadratic regression. The gray 
shaded region, nr > 0.15, is ruled out by nucleosynthesis 
constraints on the energy density of gravitational w aves with 
power law spectra IjStewart fc Brandenbergerll2008l ). We find 
that Planck by itself will not find decisive (strong) evidence 
for nr 7^ unless |tit| > 0.5 (0.43); in combination with Po- 
larBear and QUIET (II), these same levels of evidence are 
achieved for \n T \ ~ 0.36 (0.28) and 0.29 (0.25), respectively. 
These detection thresholds are ruled out for power law spec- 
tra. The constraint from nucleosynthesis should, however, be 
applied with care since the assumption of a power law tensor 
spectrum is not always appropriat^f). Meanwhile, a future 
space-based mission with the specifications of CMBPol will 
provide decisive evidence of nr 7^ for \tit\ 0.15 - just 
on the edge of the excluded region, and strong evidence for 
\ut\ 0.12. Lastly, an ideal satellite experiment will per- 
form even better: not only will it support a decisive (strong) 
confirmation of nr 7^ for \tit\ > 0.03 (0.025), but it will 
provide strong evidence for nr = 0, i.e. favor the null hy- 
pothesis Ho, if \tlt\ < 0.01. These results are summarized 
in Table 1. 



5 The Bayes factor depends on the prior range, but only weakly: 
for example, doubling the range gives |A In B01 = ln2 = 0.7. 

6 http : / / camb . info 

7 http:/ /cosmologist.info/cosmomc 

8 While power law tensor spectra are expected from slow roll in- 
flation, the primordial gravitational waves generated in scenarios 
of loop quantum cosmology are strongly non-power law: while 
typically steeply blue (ny S3 2) on CMB scales, they are scale 
invariant on smaller scales. We discuss this further in Section 5. 
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Figure 3. Bayes factor as a function of tit with r = 0.1 
for the experiments discussed in the text: Planck (black solid), 
Planck+PolarBear (magenta dotted), Planck+QUIET (II) (pur- 
ple dash dot), CMBPol (blue dashed), and ideal satellite (red 
long dash dot). The gray region is excluded for power laws by 
nucleosynthesis constraints. 



Table 1. Fiducial nj- for which the different experiments will sup- 
port strong (| In i?oi| = 2.5) and decisive (| In Soil = 5) evidence 
in favor of tlt 7^ (lnBoi < 0) and = (lnBoi > 0). 



\n T \ 



Experiment 



lnSoi 



-2.5 



\n T \ 
lnBoi = 



|n T | 
In £01 = 2.5 



Planck 
+PolarBear 
+QUIET (II) 

CMBPol 

Ideal 



0.43 
0.28 
0.23 
0.12 
0.025 



> 0.5 
0.36 
0.29 
0.15 
0.03 



0.01 



4.2 Dependence on r 

The previous findings apply to Hi with fiducial r = 0.1. 
However, the error on nr, which largely determines the 
Bayes factor, depends on the base value of r. We now ex- 
amine whether and how our conclusions change when r is 
also allowed to vary. The preferred approach to this prob- 
lem would be apply the same MCMC analysis on data sets 
generated from a range of fiducial r, in addition to nr. This 
method, however, is time consuming and more efficient ap- 
proaches exist. The previous results from the MCMC anal- 
ysis show that the posterior distributions of r(fc«) and nr 
are nearly Gaussian and not too strongly correlated, and so 
a Fisher matrix ana lysis should provide reliable constraints 
l|Perotto et alj|2006h . 

A forecast of parameter constraints can be obtained rel- 
atively easily by Taylor expanding the log-likelihood func- 
tion, ln£(0|d), about the best-fit parameter values, 0, and 



examining the 2 
d 2 hiC 



-order coefficient, 



dOidOi 



9=8 



(21) 



The Fisher information matrix, Fij, encodes parameter cor- 
relations and measures the steepness of the likelihood func- 
tion in the direction of each parameter 6i. The minimum 
precision with which para meter 8j can be measured is set 
by the Cramer-Rao bound jTegmarklll997l ). 
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(22) 



We again consider only the parameters r and w r, and only 
inclu de B-mode data in the likelihood function l|Bond et al.l 

mi), 



21n£ = ^2(2£ + l)/ sky < h> 

t 

C. 



BB + N BB 
C BB 



VBB 



+ 



C BB + N BB 



(23) 



where Cf B and N e BB are defined in Eqs. JTBJ and (|20|) . 
respectively, and C^ B is the theoretical spectrum. Using 
Eq. (23) in Eq. gives 



We now extend our MCMC analysis of CMBPol and the 
ideal satellite experiment to examine how these projections 
vary with fiducial r (Planck and the ground-based exper- 
iments fail to provide strong evidence for tit 7^ consis- 
tent with nucleosynthesis constraints even in the optimistic 
case.) We compute the Fisher matrix Eq. (|24p for each ex- 
periment on a 200x200 grid in the (riT,logr) plane with 
r £ [10" 6 ,0.4] and riT £ [—0.5,0.5]. Having obtained ay 
and Ojirp from Eq. (|22p . we obtain the Bayes factor via Eq. 
(|4)| under the assumption that the marginalized posterior 
distributions of r and tit are Gaussian. We present results 
for CMBPol and the ideal satellite experiment in Figure 4. 
The CMBPol thresholds for strong and decisive evidence 
are indicated in violet and blue, respectively, and the same 
thresholds for the ideal experiment are given in orange and 
red, respectively. For power law spectra, CMBPol will find 
strong evidence for Hi if r ^ 0.06, and decisive evidence 
if r 0.1 (in this case confirming the optimistic results of 
the MCMC analysis, Figure 3.) The ideal experiment has 
a much better chance of detecting nr 7^ 0: r ^ 10 -3 and 
r ^ 10 -4 for decisive and strong evidence, respectively. The 
ideal satellite will also find strong evidence in favor of the 



null hypothesis, Ho, for r > 2.5 x 10 
orange dotted contour in Figure 4. 



and tit within the 



5 DISCUSSION 

We have obtained projections for the detectability of tit 7^ 
with current and future satellite and ground-based CMB ex- 
periments. Current and proposed missions will struggle to 
detect a tensor tilt with statistical significance; however, an 
ideal satellite experiment will perform well. This indicates 
that discriminating between theories of early Universe struc- 
ture formation with primordial B-modes is in principle pos- 
sible using space-based platforms. We now examine the im- 
plications of our results for distinguishing between several 
competing theories for an optimistic detection of r = 0.1. 
The following results are summarized in Table 2. 

Inflation: General inflation predicts a tensor amplitude 
Ph oc H 2 at lowest order in slow roll. The dominant en- 
ergy condition demands that H < in a Friedmann Uni- 
verse, leading to a red tilted tensor spectrum, tit < 0. The 
detection thresholds for the different experiments are pre- 
sented in Table 1: while large tilts are needed for Planck 
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Figure 4. Baycs factor as a function of r and ut for CMBPol and 
the ideal CMB experiment. For CMBPol, thresholds for strong 
(lnBoi = —2.5) and decisive (lnSoi = —5) evidence in favor of 
Hi are indicated by the top two contours: violet dashed and blue 
solid curves, respectively. For the ideal mission, these same thresh- 
olds are given by the lower two contours: orange dashed and red 
solid curves, respectively, and the central closed orange dashed 
contour indicates the threshold for strong evidence (lnBoi = 2.5) 
favoring the null hypothesis, Ho- The interpretation is that in- 
side this contour rty = is strongly favored over nj- 7^ 0. The 
gray region is excluded by nucleosynthesis constraints, and the 
vertically-oriented solid black line gives the prediction of single 
field slow roll inflation. 



and ground-based experiments to make a decisive detection 
in the most optimistic case, these missions are nonetheless 
capable of providing evidence in support of inflation. Future 
missions will improve on these capabilities, but will still re- 
quire relatively large tilts in the most pessimistic outcome of 
small r. Single field inflation predicts a consistency relation, 
r = — 8nr, shown as a dashed curve in Figures 4 and 5. Our 
results indicate that current probes and proposed missions 
like CMBPol will fail to provide strong evidence in support 
of single field inflation. If single field inflation is the true 
source of the B-mode, an ideal experiment will provide de- 
cisive supporting evidence if nr < —0.02, corresponding to 
r > 0.15 (cf. Figure 4). 

String Gases: During the quasi-static Hagedorn phase 
of string gas cosmology, the Hubble radius shrinks so that 
fluctuation modes come to exist on cosmological scales 
ijNaveri et alJ I2OO6T ) . The anisotropic pressure of the mat- 
ter fluctuations generates large scale tensor perturbations; 
this pressure is smaller deep in the Hagedorn phase than 
it is during the subsequent radiation dominated expansion 
|Brandenberger et al.1120071 ). The larger scale tensor modes, 
which exit the horizon earlier, will therefore have smaller 
amplitudes than those that exit later, and the resulting ten- 
sor spectrum is a power law with a blue tilt. Because it is 
a power law, nr is bound by nucleosynthesis constraints, 
giving a small window within which rir is large enough to 
be detected, but smaller than the nucleosynthesis cutoff. 
CMBPol and the ideal satellite will provide strong observa- 
tional support if 0.12 < n T < 0.15 and 0.025 < n T < 0.15, 
respectively. 

Loop Quantum Cosmology: Motivated by loop quantum 
gravity, loop quantum cosmology (LQC) is characterized by 
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Table 2. Conditions under which the different experiments will be able to provide strong supporting evidence for different models of early 
Universe structure formation: single field inflation (SFI), string gases (SG), loop quantum cosmology (LQC), and a nonsingular matter 
bounce (MB). For LQC, we assume that nj- = 2, large enough to be decisively detected by all experiments. For the other experiments, 
supporting evidence is obtained for the fiducial values of given. These results are for r = 0.1 for all models except SFI, for which 
strong supporting evidence can only be found for r > 0.15. See the text for details. 



Experiment SFI SG LQC MB 

Planck No No Yes No 

+PolarBear No No Yes No 

+QUIET (II) No No Yes No 

CMBPol No 0.12 < n T < 0.15 Yes No 

Ideal n T < -0.02 0.025 < n T < 0.15 Yes Yes 



a quantum bounce joining a prior contracting phase to the 
expanding phase of big bang cosmology. It was discovered 
early that inverse volume corrections lead to a period of 
super inflation (w < —1) in the early Universe l|Boiowaldl 
2002), generating a power law spectrum of tensor mode s 
with an unacceptably blue tilt (|Calcagni fc Hoss ain 2009). 
More recent analyses have studied the full evolution of the 
perturbations from their generation during the contracting 
phase, across the bounce, and through a period of slow roll 
inflation driven by a m assive scalar field |Miclczarck1 l2010l ; 
iMielczarek et al.ll2010l ). Due to the shrinking Hubble radius, 
modes initiated during the contraction are strongly blue, 
Ph ~ k 2 . After the bounce, those modes that are outside 
the horizon are frozen and the blue spectrum is retained 
on these scales; modes within the horizon after the bounce 
evolve during the subsequent period of inflation leading to 
a nearly flat tensor spectrum on smaller scales. The result- 
ing spectrum is effectively a power law with a steep blue 
spectrum, tit ~ 2, on CMB scales, but the non-power law 
nature of the spectrum allows it to evade the nucleosynthesis 
constraints. This tensor tilt is large enough to be decisively 
detected by all the experiments considered in this analysis. 

Matter Bounce: String cosmology can also accommo- 
date a bouncing Universe. The early implementations, Pre- 
Big Bang cosmology and the ekpyrotic scenario, included 
a singular bounce connecting the contracting and expand- 
ing phases. Due to the dynamics of the contracting phase, 
the resulting tensor perturbations in the expanding phase 
are unacceptably blue in both models, with the requirement 
that the tensor spectrum be suppressed. Initially introduced 
as a toy model exhibiting a nonsingular transition, a generic 
bounce in which the Universe passes from matter to radia- 
tion domination prior to the boun ce is capable of producing 
a sca le invariant tensor spectrum (Fin elli fe Brandenberger] 
2008). Perturbations generated during a collapsing dust 
dominated Friedmann Universe ar e identical to those gen- 
erated during de Sitter expansion l|Wandslll999l ). since the 
amplitudes of all modes - both sub- and super-horizon - 
grow at the same rate during dust dominated contraction. 
This model is our null hypothesis, Ho, and it is phenomeno- 
logically distinct from all of the other models. If tit = 0, an 
ideal satellite will provide strong evidence for it. 

Although we have not considered them in this work, 
it is expected that future space-based laser interferometers 
like Big Bang Observer and Japan's Deci-hertz Interferom- 
eter Gravitational Wave Observatory (DECIGO) will mea- 
sure tit at least as precisely as an ideal CMB sa tellite l|Setd 
l200rj ; lKudoh~ et al . 2006; Zhao fc Baskararl2009h and should 



perform well under Bayesian model selection; we leave this 
analysis for future work. In the meantime, we can conclude 
that the awesome prospect of using the primordial B-mode 
as a window into the origin of structure formation will be 
within the grasp of next-generation space probes. 
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